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Abstract. Gravitational radiation that propagates through an inhomogeneous mass distribution is subject to random gravita- 
tional tensing, or scattering, causing variations in the wave amplitude and temporal smearing of the signal. A statistical theory 
is constructed to treat these effects. The statistical properties of the wave amplitude variations are a direct probe of the power 
spectrum of the mass distribution through which the waves propagate. Scattering temporally smears any intensity variations 
intrinsic to a source emitting gravitational radiation, rendering variability on time scales shorter than the temporal smearing 
time scale unobservable, and potentially making the radiation much harder to detect. Gravitational radiation must propagate out 
[ through the mass distribution of its host galaxy before it can be detected at the Earth. Plausible models for the distribution of 

\^ ■ matter in an L, host galaxy suggest that the temporal smearing time scale is at least several milliseconds due to the gas content 

' alone, and may be as large as a second if dark matter also scatters the radiation. The smearing time due to scattering by any 

, galaxy interposed along the line of sight is a factor ~ 10' times larger. Gravitational scattering is an excellent probe of matter 

^'^^ ■ on parsec and sub-parsec scales, and has the potential to elucidate the nature of dark matter. 
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1. Introduction 



^ Current and planned gravitational wave detectors have the potential to identify sources of gravitational radiation out to extra- 
• • , galactic distances and, for some types of sources, out to cosmological distances (Cutler & Thorne 2002) over a wavelength range 
_ ^ spanning more than 8 orders of magnitude. Instruments such as the Laser Interferometer Gravitational Wave Observatory (LIGO) 
, are sensitive to radiation in the range 10 - 10"* Hz, while the planned Laser Interferometer Space Antenna (LISA) experiment will 
' be sensitive to radiation from frequencies 10 ' - 10 Hz. Other methods of detection, such as those due to pulsar timing observa- 



tions, are sensitive to nano-Hertz gravitational radiation (Sazhin 1978, Detweiler 1979). As the detection of gravitational radiation 
requires extremely sensitive apparatus, consideration has been devoted as to whether gravitational lensing by intervening massive 
objects might focus the radiation and enhance its detection likelihood (Ruffa 1999). 

These lensing calculations (Ohanian 1974, Bliokh & Minakov 1975, Bontz & Haugan 1981, Thorne 1983, Deguchi & Watson 
1986a) are performed in the regime of physical optics. Diffractive effects are important because the wavelength of gravitational 
radiation is large compared to the Schwarzschild radius of any lensing object. For a wave of frequency v, diffractive effects are 
important for lens masses :$ lO^v"^ (Takahashi & Nakamura 2003). Thus lensing by almost any intervening single object 
must be treated using physical optics. Moreover, even for lens masses in which geometric optics is valid, it is preferable to work 
in the framework of physical optics because it allows one to compute certain quantities related to the wave amplitude that are in 
principle measurable with gravitational wave detectors but are not calculable under the framework of geometric optics. 

Foregoing lensing treatments have largely been limited to the effects of a single massive object on the wave amplitude. The 
approach has hitherto been from a deterministic rather than a statistical point of view, in that the mass profile of the lensing object 
is well specified. However, a statistical approach is more appropriate under many circumstances. 

The argument for a statistical approach stems from the large number of objects that are likely to contribute to the wave 
amplitude in many astrophysically interesting situations. There are two senses in which a statistical approach is applicable to the 
gravitational lensing of gravitational radiation. One is when a large number of lensing objects can drift in front of the line of sight 
to a source, so that an average over time would involve a sufficiently large number of objects to merit an average over an ensemble 
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of lensing objects. A second, more restrictive, sense applies when a large number of objects contributes to the observed wave 
amplitude at any instant. This latter sense is likely to apply in practice. The wave amplitude of a scattered signal is determined by 
the phase fluctuations caused by lensing objects in a region at least as large as the first Fresnel zone (see, e.g., Narayan 1992). For 
gravitational lensing occurring at distance D from the observer, the Fresnel radius rp ~ -\JcD/2jtv, is of order a parsec for lensing 
of V ~ 1 Hz radiation at cosmological distances. Thus the area defined by the first Fresnel zone encompasses contributions from 
a large number of objects. Consider, for example, the lensing influence of stars on the gravitational radiation passing through the 
spiral arm of a galaxy. Even the volume of a line of sight orthogonal to the disk axis of a spiral galaxy with radius rp ~ 1 pc 
would typically encompass many thousands of stars. The argument for a statistical approach is strengthened when one considers 
the possible contribution of gas and of dark matter, both of which may be inhomogeneous on sub-parsec scales. 

Scattering gives rise to a number of important effects that influence the detectability of signal. The two that motivate the 
present treatment are (i) the focusing and defocusing of radiation by phase fluctuations, leading to temporal variability in the 
observed wave amplitude, and (ii) the temporal smearing of intrinsic intensity' fluctuations in the source. 

Temporal fluctuations in the wave field occur whenever there is movement of the lensing system relative to the line of sight to 
the source of gravitational radiation. Fluctuations in the gravitational potential drifting transverse to the line of sight induce both 
phase and amplitude fluctuations in the gravitational radiation, leading to temporal decorrelation of the wave field on sufficiently 
long time scales. This effect is important if the decorrelation time scale associated with gravitational scattering is short compared 
to the time scale on which the wave field would oscillate in the absence of scattering. In this paper we concentrate on fluctuations 
in the wave amplitude itself, while in a second paper we consider fluctuations in the intensity. 

Temporal smearing is due to the fact that the radiation reaching an observer can arrive from a number of different angles. 
This effect is significant when the phase fluctuations induced by the lensing are sufficiently large that radiation from outside the 
Fresnel radius is scattered toward the observer. The radiation then arrives from a number of different regions in the scattering 
region, in an effect known as multipath propagation. There is a range of propagation times associated with the range of angles 
over which the radiation arrives. The effect is important because it potentially smears out the intrinsic intensity variations of a 
source. For instance, time variations of the source would be undetectable if the intrinsic variability time scale were significantly 
shorter than the scatter broadening time. 

The scattering of gravitational radiation potentially constitutes an exceptional probe of the distribution of matter out to the 
distances to which sources can be detected. In this context it is useful to compare gravitational scattering to the interstellar scat- 
tering (ISS) of radio wavelength radiation by the diffuse ionized component of our Galaxy's interstellar medium. Both scattering 
phenomena require a statistical approach and occur in a regime in which physical optics is important. ISS yields information on 
the power spectrum of turbulent fluctuations in the ionized plasma distribution of our own Galaxy over a range of scales spanning 
some five orders of magnitude in wavenumber (Armstrong, Rickett & Spangler 1995). This information is gleaned primarily from 
pulsars which are effectively point sources, and thus constitute excellent probes of fluctuations in the interstellar medium. Other 
radio sources are less subject to ISS, and their larger angular diameters severely truncate the range of wavenumbers over which 
the power spectrum can be measured. Gravitational wave sources are the analogues of pulsars in the regime of gravitational 
scattering. They emit coherently and possess small angular diameters so the source structure does not influence the scattering 
characteristics, unlike most sources of electromagnetic radiation that are subject to gravitational lensing. Thus the scattering of 
gravitational radiation is a fine-scale probe of the distribution of dark matter in the local universe. 

The close relationship between gravitational and interstellar scattering implies that many phenomena identified in this paper 
are analogous to effects previously identified in the context of interstellar scintillation. For example, the temporal broadening 
of gravitational radiation is analogous to the temporal broadening of pulsar radiation, and it poses similar limitations on source 
detectability (see the review by Rickett 1977 and references therein). 

The outline of the paper is as follows. In Sect. 2 we review the effect of an arbitrary mass distribution on the amplitude of a 
gravitational wave. In Sect. 3 we compute the temporal variability of the wave amplitude by considering the power spectrum of 
the wave amplitude fluctuations. The effect of temporal broadening on the scattered signal is considered in Sect. 4. Estimates of 
the power spectrum of mass fluctuations due to the various constituents of a galaxy are presented in Sect. 5. In Sect. 6 we estimate 
the magnitude of scattering effects and discuss their relevance to the observability of gravitational waves. The conclusions are 
presented in Sect. 7. 

2. Propagation of the wave amplitude 

We briefly outline the equations governing the propagation of gravitational radiation in the gravitational potential of a distribution 
of lensing objects. One writes the spacetime metric in the form 

ds^ = -(1 - 2U(r))cW + 2U(r))dr^ = gl^^^dx^'dx^ (1) 

where U{r) <K 1 is the gravitational potential of the lensing objects. Consider a linear perturbation hf,y in the background metric 
tensor, g^^y-. g^y - g^^y + hf^y. If the wavelength of the gravitational radiation is much smaller than the typical radius of curvature 

' Throughout the text we refer to the quantity defined by (p{v) in Eq. J4j as the (complex) wave amplitude at frequency v, while the quantity 
I0(i')p(f) is the intensity of the radiation at frequency v and time f. Thus intensity variations refer to changes in the quantity |0(v)|(r) with time. 
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of the background metric, then one has 2R^^\ /i"^ = 0, where R^^\ is the background Riemann tensor, and one identifies the 
gravitational wave as 

hfiv - <P ^iiv^ (2) 
where e^v is the polarization tensor of the gravitational wave (ej^ = 0, e^^ys''^ = 2) and = 0(f, r) is the scalar wave amplitude. 
This scalar wave propagates according to ^f-^^g'^^^^^dycf)) - 0, which, when combined with equation ([Q, takes the following 
simple form 

{■^^ + a?)^ ^ Aa?U^, (3) 

where o) - 2nv and ^(v, r) is the temporal Fourier transform of the scalar wave amplitude. 

We consider the solution of Eq. (|3} with reference to the lensing geometry shown in Fig. 1 . A point source of unit intensity 
is located an angular diameter distance Ds from the observer and distance Dls from the lens plane. The lens plane is located 
an angular diameter distance Dl from the observer's plane, and the observer is at location X' on this plane. The gravitational 
perturbations are assumed to be located on a thin lensing plane. One might be concerned that the thin lensing approximation 
used here, while applicable to the lensing of electromagnetic radiation, is not applicable to gravitational radiation because of its 
short wavelength. The depths of most lensing objects are small compared to the wavelength of gravitational radiation. However, 
the thin lens approximation is still valid because it depends only on the weakness of the wave amplitude and not on the short 
wavelength assumption (Thorne 1983). 

The wave amplitude has the solution (Schneider, Ehlers & Falco 1992) 



0(v,X') = 



Dsv(l +Zl) 



-m/2 



cDlDls 

where the time delay is, apart from a constant, given by 



J"<i^xexp[27r/vfj(x,X')], (4) 



..(x,X') = ^^iili^ 
cD,D,. 



(5) 



The term Dis IDs in front of the observer's plane co-ordinate X' is a correction factor that accounts for the spherical nature 
of the wavefront in the mapping of lens plane co-ordinates to co-ordinates on the observer's plane (see Goodman & Narayan 
1989). When the source is located sufliciently far behind the lens plane the incident wavefront may be regarded as planar and this 
correction factor tends to unity. We introduce the scaled co-ordinate X - X'Dis IDs to retain the symmetry between lens and 
observer plane co-ordinates below. 

The phase delay associated with the surface mass distribution, S(x), is 

T« = 1^ f,.,.„,.)ln(t:^), (6) 

Dsc^ J \ xo I 

where xi) is an arbitrary normalizing constant that is of no interest because it only introduces an arbitrary constant phase delay; 
we henceforth write xq = 1 to set this phase delay to zero. 

It is convenient to recast the Fresnel-Kirchoff integral @ in the form 



-ot/2 r \ i r /lx-x'l\ 

4>(v, X) = d^x exp — T (x - Xf + iK d^x'-L(x') In ■ ■ 

2nr^ J [2rp J \ xq ) 



where K = -Sn^^^^, (7) 



and where we define the curved-spacetime generalization of the Fresnel scale in Euclidean space as 
2 _ DlDlsA 



2nDs(l +zl) 

and the phase delay due to the gravitational potential is identified expUcitly as 



(8) 



iA(x) J c/V 2(x')ln^- 



(9) 

xo 

The approach adopted in this paper involves regarding the gravitational phase delay, tfr, as a random variable. We define r^sf as the 
length scale on the lens plane over which the root-mean-square phase delay changes by one radian. The magnitude of eff'ects due 
to gravitational scattering are quantized by comparing the length scale of gravitational phase fluctuations, rdifr, to the length scale 
of geometric phase fluctuations, characterized by the Fresnel scale, rp. Scattering effects are important when rjiff is comparable 
to or smaller than rp. To see this, consider the generic properties of the wave field as determined by Eq. Q. The dominant 
contributions to the wave field come from regions where the argument of the exponential in Eq. varies slowly as a function of 
X. For small gravitational phase fluctuations on the length scale rp (i.e. rjiff » rp), only the region around |x - X| < rp contributes 
to the wave field, and mild phase curvature on the scale of rp gives rise to mild focusing and defocusing of the gravitational 
wave. For very large phase fluctuations (i.e. r^ff <sc rp) many separate regions on the lensing plane contribute to the observed 
wave amplitude. This is because there are many regions on the lens plane for which the exponential term in Q varies slowly as 
a function of x, since the large phase fluctuations can offset the contribution of the (x - X)^/2rp term. 
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lens plane 



Fig. 1. The lensing geometry. Co-ordinates on the lens plane are denoted with x and those on the observer's plane by X. 

2.1. The effect of finite source size 

As sources of gravitational radiation are compact and spatially coherent it suffices to treat the radiation from these systems 
as point-like (Mandzhos 1981, Ohanian 1983, Schneider & Schmid-Burgk 1985, Deguchi & Watson 1986b, Peterson & Falk 
1991). Thus equation suffices as a practical description of the wave amplitude observed from any source of gravitational 
radiation. Equation also describes the gravitational lensing of electromagnetic radiation from compact objects, but it needs to 
be modified once the source angular diameter becomes sufficiently large. In practice, most sources of electromagnetic radiation 
are not spatially coherent and are sufficiently large that effects due to their finite angular diameter are important. 

The effect of finite source size on the temporal variability of gravitationally scattered radiation is not discussed in the present 
paper, but the results in this paper are relevant to the lensing of electromagnetic radiation when the source size does not exceed the 
characteristic angle through which scattering deflects the radiation. Stated another way, source size effects are important when the 
source angular diameter exceeds the angular scale over which the wave field of a point source would vary due to inhomogeneities 
in the scattering medium (e.g. Little & Hewish 1966, Salpeter 1967). To see this, one can attribute wave field decoiTelation length 
scales to both the effects of source size and of scattering. Define the length scale on the observer's plane, Xsrc, over which structure 
in the source causes the wave field to decorrelate (i.e. (0(Xsrc + x')0*(x')) declines on the length scale Xsk). For an incoherent 
source of angular size flsrc, this length scale is X^rc = ^/^nOsrc- One can define a similar length scale Xscat - ^/27r0scat over which 
the wave field from a point source would decorrelate due to scattering. Source size effects are therefore unimportant when the 
source decorrelation scale exceeds the scattering decorrelation scale, ^ Xscat or, equivalently, when the angular scale of the 
scattering pattern from a point source exceeds the source angular diameter, 0scat ^ ^src- 

2.2. Statistics of the mass fluctuations 

In the following discussion we are interested in the statistical properties of the wave amplitude as influenced by gravitational 
lensing. We explicitly ignore (non-stochastic) large-scale phase gradients due to the overall mass distribution in a lensing system. 
Large scale phase gradients contribute to macrolensing, but are not relevant to the effects under consideration here. 

Since the statistics of wave amplitude fluctuations depend on the statistical properties of the phase, i^(x), we must characterize 
the fluctuations in the mass surface density, E(x), that drive the phase perturbations on the lensing plane. The two most useful 
quantities are the mass surface density autocoiTelation function 



Cdr) = <A2(r' + r)AS(r')) = <[£(r' + r) - £][E(r') - £]), 



(10) 



and its associated quantity, the mass surface density structure function. 



D^(r) = <[2(r' + r) - E(r')]') = 2[C^(0) - C^(r)], 



(11) 



which describes the mean square difference in the mass surface density between two points separated by a vector r on the lensing 
plane. The autocorrelation function is related directly to the power spectrum of mass surface density fluctuations '^^{qx, Qy) as 
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follows: 



CE(r) 



/^C>.(q).'-. (12) 



These average quantities are written as functions of separation r only, because it is assumed that the statistical properties of E do 
not depend on the position where they are measured, which is to say that the statistical properties of 2 are wide-sense stationary 
(see, e.g., Mandel & Wolf 1995). Although this assumption is commonly used in scattering theory, it is obviously flawed to some 
degree because the statistical properties of the lensing matter are expected to vary as, for instance, one moves from a line of sight 
intersecting the centre of a galaxy towards one that intersects its rim. However, it is an excellent approximation in practice since 
we are concerned with the statistical properties of the matter distribution mainly on scales less than ~ 10 pc. 

It is often convenient to relate correlations in the gravitational phase i^(x) to the power spectrum of gravitational mass fluctu- 
ations. The autocorrelation in phase fluctuations ^i'(x) and ^i'(x + r) is derived using the identity, 

^ t/2aln|a|e'i" = -2nq-^, (13) 
which yields the following relation between phase autocorrelation and the mass power spectrum, 

<AiA(x)Af (X + r)) ^K^ ^ d^^d^On III In ICI C^^ - C-r) ^ jd\ e'l-" q-' <l.s(q). (14) 
The structure function of phase fluctuations D^(r) is related to the mass power spectrum according to 

D^{r) ^2K^ J d\ [l - e'"--] q'^ Os(q). (15) 
These relations simplify the calculations performed below and allow us to express the results in physically useful quantities. 



3. Wave amplitude fluctuations 

Given the exact mass distribution on the lens plane, eq. Q can be evaluated directly to find the wave amplitude on the observer's 
plane. The mass distribution on the lensing plane is usually unknown, but it is possible to express the statistical moments of the 
wave amplitude in terms of statistical properties of the mass distribution. In the following subsection we evaluate the variability 
due to gravitational scattering by computing the covariance of wave amplitude fluctuations across the observer's plane. 

3.1. Mean visibility 

We wish to compute the time scale on which the wave amplitude of a source lensed by random mass inhomogeneities fluctuates. 
The temporal fluctuations are characterized by the covariance between the wave amplitude received at a frequency v and time f 
with those detected at the same frequency at some later time t + At: {4>{v; t)4>*{v; t + Af)). The angular brackets here refer to an 
average over all possible statistical ensembles of lensing objects which, here, is equivalent to an average over all times t. 

The mass fluctuations on the lensing plane are assumed to be fixed, and the lensing plane advects past an observer at some 
effective velocity Veif . This is known as the frozen screen approximation. A stationary observer who samples the wave field at 
position X at times fi, f2, . . . would observe exactly the same wave field as an observer who moves across the observing plane, 
measuring the wave fields at all the positions X + Vefffi, X + Vgfffa^ ... on the observing plane. An average over time is equivalent 
to an average over all positions X on the observer's plane. To be specific, the temporal covariance of the wave field is equivalent 
to the covariance measured between two separated receivers located at positions x and x -H X = x -H Vetf Af . The specific expression 
for the effective velocity, in terms of the velocities of the source, lensing plane and observer is (see Appendix|B} 

Veff = Vscreen - VEaith ( 1 " | " Vsrc ( ^ | ■ (16) 



To compute the wave field covariance we consider the wave amplitude measured by two detectors a distance X apart on the 
observing plane: one receiver is positioned at -X/2 while the other is located at X/2. In light of its connection with the corre- 
sponding quantity in interferometry (see Thomson, Moran & Svenson 1986), we refer to this covariance as the mean visibility, 
or mutual coherence, of the gravitational wave field. Using the the Fresnel integral derived in eq. (0, the mutual coherence is 

y(y; X) = <0(v; -X/2)^*(y; X/2)) = —L^ f d'-xd^x' exp [(x + X/2f - (x' - X/lf] 

+iK j flf24S(a[ln||-x|-ln||-x'|]|. (17) 
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It is useful to make the change of variable r = x - x' and s = (x + x')/2, so the visibility simplifies to 

iS-(r + X) 



y(y;X) = -i— r dhd'ssxp 



+ iK d"-^ In 1^1 + r/2 + s) - 2(| - r/2 + s)] 



(18) 



We assume that the mass-induced phase fluctuations are wide-sense stationary, so that the difference of the mass surface density 
+ r/2 + s) - - r/2 + s) depends only on the separation r. The integral over s then yields a factor (2nr^)^S^{r + X) and 
after performing the remaining integral over r, the visibility reduces to 



V(v;X) 



|exp iK J t/^mi) (In II -I- X/2| - In II - X/2) 
The average over the mass fluctuations is evaluated using the result 
<exp[ia]) = e'<"> 



1 + i{6a) - i<fc2) _ i.(5«3^ ^ 



(19) 



(20) 



for a random variable a with average {a} and fluctuations 6a. This approximation makes no specific assumptions about the 
distribution of the random variable da. When the fluctuations in the gravitational phase, i/', are small, it is sufficient to consider 
only the first three terms in this expansion: 



<y(y;X)> * I + iK{Y} - ^K\y\ 



where Y 



I In III [m-X/2)-S(| + X/2)]. 



(21) 



The imaginary component of the contribution is zero because {Y} = 0; for mass fluctuations which are wide-sense stationary (i.e. 
their statistical properties are independent of position) the average (Z(| - X/2)) is equal to (Z(| + X/2)). The final term involving 
terms 0{K^) simplifies to 



X 



iai>o,|/3|>o 



d^OLd^P \n\a.\ ln|/3| 



D^{a -(B-X) D^{a -j3-X) 



£>!(« - /3) 



2 2 

which, re-expressed in terms of the power spectrum of mass surface density fluctuations using equation il4\ . yields 
<y(y, X)}^l-K^ J d\ q-' 0£(q) [1 - cos(q ■ X)] . 



(22) 



(23) 



One obtains a clearer result on the length scale over which the scattered wave field varies by examining the power spectrum 
of wave amplitude fluctuations across the observer's plane, which is just the Fourier transform of the visibility: 



J d^Xe-"i''{V(v;X)} = KY^(t>^(l) + (^nf^^l) 1"^'/ dW q'-%^(q') 



(24) 



Apart from the unimportant zero-frequency (q - 0) contribution, the spatial Fourier transform of the mean visibility yields a 
direct measurement of the power spectrum of mass density fluctuations weighted by q^*, which weights the measurement to the 
largest lensing structures present in the medium. Measurements of the mean visibility potentially provide an elegant probe of the 
mass density fluctuations along the line of sight to a source of gravitational radiation. 



3.1 .1 . Strong phase fluctuations 

The foregoing result embodied in equation (I24> is only correct in the limit in which the phase fluctuations are small. Using 
equation il5\ . we can express the visibility in terms of the phase structure function: {V(v, X)) ^ 1 - D^^,(r)/2. This shows that the 
small phase approximation is only valid for sufficiently small baselines r such that the mean square difference in the phase delay 
due to random mass fluctuations across the lensing plane is less than one radian. 

It is possible to employ an alternative approximation to derive the mean visibility in a manner that is independent of the 
magnitude of the phase fluctuations. However, this requires that we make specific assumptions about the statistical properties of 
the gravitational phase, ifr. In particular, we assume that phase fluctuations are normally distributed, which allows the average 
over phase fluctuations to be performed using the result 

<exp[!iA]) = exp[!XiA> - <(5<A')/2]. (25) 

Applying this result to the phase fluctuations in ( I19> . one derives an exact result for the mean visibility under this somewhat more 
restrictive assumption: 



y(v;X) = exp 



[ £'^(X)1 




2 


= exp 



j d\q-*[l-e''i-^] cE)j;(q) 



(26) 
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Fig. 2. Temporal smearing of a time variable signal (i.e. one for which |0(y)|(f) varies with time) smooths out the variations so 
that the radiation arrives at the detector over a longer interval and the source potentially falls below the detection threshold. 

This generalization of (I23> shows that the wave amplitude on the observer's plane is only correlated over distances on which 
the mean square phase difference on the lens plane is less than one radian. In order words, the wave field is only correlated over 
a distance ~ rdiff. (Recall that raiff is defined such that Z)^(,(rdiif) =1.) The frozen screen approximation described above then 
implies that temporal fluctuations in the wave field occur on time scales as short as fvar ~ '"diff/i'eff- The covariance between wave 
amplitudes measured on time scales longer than fvai decreases exponentially quickly. 

We can exploit the similarity between scattering due to gravitational mass perturbations and interstellar scattering in inter- 
preting the physics embodied in equation ( I26l l. For phase fluctuations 0(x) in the interstellar medium, the visibility of a scattered 
point source is yiss(T, X) - e\p[-D^(X)/2], which is identical in form to equation i26\ . The visibility is related to the apparent 
source brightness distribution (i.e. the source image), I{6), according to (e.g. Thomson et al. 1986) 



Since the average visibility of a scattered source declines on a length scale X — rdiff, the average angular brightness distribution 
of the scattered source declines on an angular scale - l/kras. The implication for observations of gravitational radiation is 
that a detector observing a scattered source receives radiation from a range of angles out to an average angular-broadening scale 
6'scat = i/krdiff. 

Under certain circumstances it is appropriate to consider the variance in the visibility as a measure of the extent to which it 
deviates from its mean value. The fourth moment of the wave field, which describes the covariance of fluctuations in the visibility, 
is computed in a subsequent paper (Macquart, in prep.). 

4. Temporal broadening of the scattered signal 

Propagation effects can cause temporal smearing of gravitational radiation, leading to possible misinterpretation or even non- 
detection of bursty phenomena. Temporal broadening is important when a signal propagating through a random medium can 
potentially reach an observer along several different paths (this is known as multipath propagation). As different time delays are 
associated with the ray paths, power from an impulsive event is smeared out in a process known as temporal broadening. This can 
redistribute power in the emitted signal so that the temporally varying signal is smeared to a level below the detection threshold 
(see Fig. 2). 

The temporal smearing of gravitational radiation is analogous to the temporal smearing of radio wavelength radiation re- 
ceived from pulsars. For pulsar radiation this effect is due to scattering off electron density inhomogeneities in the interstellar 
medium, and is important for most pulsars located in the Galactic plane, as it limits the temporal resolution with which the pulse 
shape can be measured. The effect severely limits the detectability of pulsar radiation when the temporal smearing time scale is 
comparable to the pulse period. This is particularly important for the heavily scattered pulsars located near the Galactic centre, 
as the broadening time scale is so large that pulsar emission appears continuous rather than pulsed; although the emission is still 
detected, it is not time- variable. 

The treatment developed for the smearing of gravitational radiation is analogous to that derived for the temporal broadening 
of pulsar radiation due to interstellar scattering (Lee & Jokipii 1975). The main differences here are that perturbations in the 




(27) 
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gravitational potential drive the phase fluctuations rather than fluctuations in the electron density, and that the gravitational phase 
fluctuations scale like /I"', whereas those in an interstellar plasma are proportional to A. 

Before discussing eff'ects caused by inhomogeneities in the gravitational potential, we first consider effects related to the 
propagation of the wave through a homogeneous gravitational field and its subsequent detection. We explicitly decompose the 
wave function into its temporal Fourier components, 0(w) as follows: 



cf>(t) 



1 r 

27r J-o< 



cf>{tj)e 



' do). 



(28) 



The first effect that limits the measurement of the wave function is the spectral response, or bandpass, of the detector, which we 
idealize as 



= exp 



-(OJ - (jJqY 



2A2 



(29) 



(30) 



where wo is the central frequency and A is the bandwidth. The detected signal is then 

4>{co)fB(co)e-""daj. 

Consider the response to radiation with a sharp pulse at time f = propagating through free space. Then (p(t) - 6(z/c - t), 
4>ico) - e'*-" and the measured wave function is 



I0det(f)l = 



exp 



(f - z/cfA^ 



(31) 



Now we introduce effects related to propagation through a homogeneous gravitational field. Consider first a homogeneous 
gravitational field which gives rise to a constant phase delay (if/). One can relate this to a decrease in the group velocity of the 
radiation with respect to c. A phase delay iff gives rise to a time delay ifr/cj and, for propagation through a distance z the effective 
group velocity of the radiation is - (dk/dcoy^ - c — zoj/i//, and J31> becomes 



I0det(f)l = 



exp 



(t-z/v.^A^ 



(32) 



Phase perturbations due to gravitational inhomogeneities are not dispersive because the group velocity is independent of fre- 
quency. As a result, there is no additional decorrelation due to any dependence of the effective refractive index of the lensing 
medium on frequency. 

We now investigate the effect of an inhomogeneous mass distribution on the temporal broadening of the signal. Using the 
decomposition given by equation ( I28t . the time-varying intensity may be written as 

(33) 



1 r°° 

{2ny J_co 



where the wavenumbers k and k' are functions of oj and cj' respectively. The terms inside the integral depend, in general, 
on both the frequency a> and the frequency difference oj - co'. However, if the bandwidth of the receiver is small, the quantity 
{(pik, z, r)^*{k', z, r)) depends strongly only on the frequency difference Aw = u-u' , and its dependence on the central frequency 
h) (or w') can be neglected. One then writes Aw - (jj - u' - {u + co')l2 and performs the integral over W, leaving 



/det(0 



(2 



c/Aw {^(aj)^*((L) + Aoj)} 



This is recognized as a convolution with respect to /, so that 

/det(f) = ^l(f)*^2(f), 

with 



and 

Piit) 



o, 



dAk{^(k)^*(k + Ak)}e 



icAk(t-z/Vg) 



(34) 



(35) 



(36) 



(37) 



In practice the convolution described in ( I35t means that either Pi or P2 is primarily responsible for the temporal broadening, 
depending on which of the two is the broadest function of t. The term Pi represents the temporal intensity variation observed 
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from a pulsed signal after propagating through a homogeneous gravitational potential. In the present case, when the detector 
bandpass is large, this term takes the limiting form 6{t - z/g), which represents the arrival of a sharp pulse at time f = z/vg. More 
generally, when the detector bandpass is large we can write the temporal intensity variations as the convolution of the intensity 
variations intrinsic to the source, Isidt), with the term P2(t)- 



/det(0 = hvc(t) * ^2(0- 



(38) 



The term f 2 embodies the effect of mass inhomogeneities on the temporal broadening. This term is the Fourier transform of the 
visibility between the wave field measured at some wavenumber k\ - k with the same field measured at wavenumber k2-k + A^, 
and with X = 0. Using equation (0), this mean visibility is 



Ak 



(39) 



where we have written ki - k + Ak/2 and k2 = k - Ak/2. The quantity {4>{ki,z)4>*(k2,z)) describes the spectral decorrelation of 
the signal between wavenumbers k\ and ^2 caused by its propagation through the inhomogeneous gravitational potential between 
the source and the observer. We assume that the phase fluctuations obey gaussian statistics in order to compute the average over 
the phase fluctuations. However, in interpreting the physical origin of effects due to temporal broadening, we shall see the results 
obtained here are quite generic, and the assumption of gaussian statistics has little bearing on the nature of the effects. We also 
neglect terms of order Ak^ jk^ or higher except in front of the potentially large C^(0) term. The mean visibility evaluates to (see 
AppendixlAl 



V(k 



>^2,z) = £/^rt/2qexp 



-ir ■ q + 



iq^Akrl Ak^ 



2k^ 



k^ 



(40) 



Equation ( I4Q> is identical in form to the equivalent expression found by Lee & Jokipii (1975) for scattering broadening in the 
interstellar medium due to a thin layer of electron density fluctuations, and it is convenient to follow their treatment of the effect 
of the phase fluctuations on the profile of a temporally broadened signal. Replacing Eq. J40l i in Eq. J37l i. P2 can be represented 
as the convolution of two broadened profiles 



Piit) = Puit) * Poit), 
where 



dAk exp 



Ak^ 

-icAkit - z/vg) - —C^(0) 



and 



Poit) 



(2: 



dAk d^q d^r exp 



-icAk(t - z/vg) ■ 



iAk J , 

'■r ■ q + —1 4 ■ 



^O^(r) 



(41) 

(42) 
(43) 



We refer to Pr and Pd as the refractive and diffractive broadening terms respectively, for reasons explained below. Performing 
the integrals over Ak, one obtains a closed expression for Pr and writes Pq in terms of a Fourier transform: 



PR(t) 



ck 



Poit) = [ dhd\5 

l,J^^rexp 



4Q(0) 



and 



2kc 



exp 



DJr) 



-ir ■ q - 



DJr) 



-ir ■ q - 



where q — — 



cklt 



(44) 



(45) 



(46) 



The effect on the pulse shape due to these two terms is qualitatively different. The pulse broadening due to the refractive term is 
a gaussian symmetric about the mean pulse arrival time with width 



ck 



^/q(0). 



(47) 



The symmetry of the refractive pulse broadening profile implies that the pulse can arrive either earlier or later than the mean 
arrival time t - z/vg. One understands the origin of this effect in terms of deviations in the total gravitational potential along the 
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observer 



source 



Dls 



lens plane 



Fig. 3. Diffractive temporal broadening arises because radiation, scattered through a typical angle flscat = (^'"diff) ' , is delayed 
relative to radiation propagating directly along the Une of sight to the source. 

line of sight. These give rise to deviations in the total phase delay along the line of sight about the mean value (i^). The time 
delay due to a deviation of phase Sif/ is just di/z/ck, so the root mean square phase delay expected from variations in the total 
gravitational potential is {6if/^y^^/ck, which is exactly the temporal broadening time given above. Refractive broadening only 
affects the mean pulse profile over a long time (see §6.1), and is not relevant to the effect depicted in Fig. 2. This is because the 
broadening occurs only on a time scale over which the total phase delay along the ray path varies, which is, in practice, large. 

On the other hand, pulse broadening due to the diffractive term is asymmetric: inspection of Eq. ( I46> shows that ^^(f) is zero 
before the pulse arrival time f = z/vg. Unlike refractive broadening, diffractive broadening affects the shape of each individual 
pulse, and is directly applicable to the effect described in Fig. 2. The exact form of the diffractive broadening term has not been 
determined analytically, except for the specific case - (r/raiff)^, where one derives 



The temporal broadening time scale is estimated for other phase structure functions by inspection of equation (I46> . The diffractive 
profile begins to decline when 1 /q is equal to the scale on which the term exp[-D^(r)/2] varies. This occurs on the scale r = raiff, 
so Poit) declines on the time scale given by -sjck{t - z/vg)/rf = r^^. Thus the temporal broadening time is 



The origin of the diffractive broadening term lies in the fact that the observer receives radiation from a range of angles over the 
lensing screen. Radiation that is scattered through a larger angle takes longer to reach an observer. The typical scattering angle 
is fiscal - l/^'"diif ■ For an observer located a distance z behind the lensing plane the time delay associated with radiation arriving 
from angle l/kruff is z(ck^ r^^^y^ , which is equivalent to i49i above. The geometry associated with this time delay is illustrated 
in Fig. 3. Numerical estimates of the scatter broadening profile due to the diffractive term are given in Lee & Jokipii (1975). 

5. Models for the mass fluctuations 

Since the foregoing results depend on the power spectrum of the surface mass density fluctuations, it is appropriate to quantify the 
form and magnitude of this spectrum. We consider the power spectrum of mass fluctuations in a galaxy due to its stellar, gas and 
dark matter content. Although it might be supposed that dark matter dominates the mass content of a galaxy, the extent to which 
it contributes to gravitational scattering relative to other forms of matter depends critically on its distribution. With the nature of 
dark matter largely unknown, the power spectrum of dark matter fluctuations on scales relevant to those probed by gravitational 
scattering is ill-constrained at present. Nonetheless, it is necessary to consider the conditions under which the effects are likely 
be measurable, even if the estimates are rough. 

The results derived here are employed in the following section to estimate the magnitude of the effects caused by gravitational 
scattering. 

5.1. Scattering by a star field 

Here we consider the mass power spectrum due to a collection of stars, each of mass M and with positions ri, r2, . . . , r^. The 
mass surface density is written in the form 




(48) 



N 

2(r) = 2M/(r-r,), 

i 



(50) 
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where Mf(r) is the projected surface mass density of a star centred on the origin. We compute the mass surface density covari- 
ance: 

<S(r' + r)S(r')> = j J dh'i:(r' + r)2(r') = j d\e-''i'\t(q)f. (51) 

The power spectrum of the surface density is 

N 

|£(q)|2 = ^ M\f{qfe'^<'''-'^ (52) 

ij 

and its average is computed by separating the self (/ = j) contributions from the A^(A^ - 1) cross-term (/ j) contributions: 
\t(q)f = NM^imf + N(N - l)MV(q)l' {e"'-^"'-'^^) ■ (53) 

The average over positions r, and appearing in the last term vanishes if the masses are distributed randomly over all space. 
However, this average is in general non-zero and one has 

^^-■q(r,-r,)^ ^ i_ J d'r,d'rjp2ir,,rj)e"i-''^-'^'\ (54) 

where /?2(r,, r^) is the joint probability of finding an object at position r, while another is located at position r^. In practice this 
average is never zero because the idealization of stars distributed randomly on an infinitely extended lensing plane is unrealistic. 
The objects are always confined within a finite area, however large, so there is always some maximum separation between pairs 
Tmax = r, - Tj. Thus there is always a sufficiently small q ~ l/r^ax for which the average in J54> is equal to one. However, 
in certain situations it is convenient to separate out the large-scale distribution of matter, which may give rise to large phase 
gradients and hence gravitational macrolensing, from the local variations in mass density due to stars. In other words, it is often 
convenient to treat the distribution of stars as being locally uniform, even when the stellar distribution cannot be uniform on large 
scales. 

Now if the distribution of object positions depends strongly on the difference, Ar - r, - r^, and only weakly on the average 
position s = (r, + Vj)/2, we can approximate the distribution of object positions, as described by the joint probability p2, as being 
wide-sense stationary. We write p2iri, Vj) as a function of Ar only, which gives 

(e'q(r,-r.)^ = ^2 / d^^rd^Asp2{Ar)e''i'"' = (55) 



Placing these averages back into equation ( I5U and, using the fact that the mean surface density. (2) is (N/A)Mf(q = 0), we 
obtain the average autocorrelation of the mass surface density fluctuations 



j d\e~''^-'\f(q)Hl+crp2(q)) 



cr^M^PiO), (56) 



(27r)2 

where cr - N/A is the mean surface density of objects. The corresponding power spectrum of surface density fluctuations is 
(D2(q) = crM^imf [1 + crp2(q)] - cr^M^6(q)f(0). (57) 

This power spectrum is inserted directly into equation il5\ to determine the structure function of phase perturbations due to 
gravitational lensing by a field of stars. 



5.1 .1 . Scattering by a homogeneous stellar distribution 

We construct the phase structure function due to a population of objects with mass M. In order to construct a model 
that is analytically tractable, we idealize the radial density profile of the objects as a gaussians with radius R: p{x,y,z) = 
M (IttR^)-^^- exp[-(x2 + z^)/2R^]. The surface density of each object is then S = M (2nR^)-^ exp[-(x^ + y^)l2R\ from 
which we identify (cf. eq. \50V ) 



where r - ^Jx^ +y^, q - -^ql + q^ (58) 

When the outer scale of the stellar distribution is much larger than the scale sizes of the lensing objects, one separates the 
contribution due to the large density gradient introduced by the overall stellar distribution from the local mass density fluctuations 
caused by individual stars. If the stars are distributed homogeneously locally, then the power spectrum of local mass fluctuations 



/(r) = 



2nR^ 



exp 



2R^ 



and /(q) = exp 



2n2 



q^R 
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reduces to <l>2(q) = crM^|/(q)p, where we neglect the unimportant term proportional to 6(q) in eq. ( I57t . The phase structure 
function, D^, is constructed by substituting the power spectrum of mass surface density fluctuations in Eq. il5\ : 



D^{r) = 'XnK^a-Nf J dqq'^ [1 - Jaiqr)] exp [-q^R^\ ■ 



(59) 



The q^^ weighting term of the integrand causes problems with the convergence of the structure function as ^ — » 0, because the 
rest of the integrand does not increase sufficiently quickly at small q. It is convenient to separate part of the integrand, 1 - iaiqr), 
into two terms (1 - Joiqr) - q^r^/4) and q^p- jA. With this separation, one has 



h 



r -3 

= 1 dqq 

Jo 

4 9'-o I , 



1 - J^iqr) — 



(2-r)r2 



dqq 



HmEi i-q^R^) lim(r + logC^^^^^) , 



(60) 
(61) 

(62) 



The integral /i converges because the term (1 - Jo(qr) - q^r^/4) contains only terms Oiq"^) and higher. The integral I2 diverges 
because it contains terms 0{q^^). However, as can be seen in eq. i6H . the divergence is only logarithmic in q, so we approximate 
the contribution due to I2 as multiplied by a term of order unity. Physically, this is justified because there must always be some 
large scale at which the power spectrum of mass density fluctuations cuts off' because the distribution of stars is finite. In practice 
the stars are only distributed homogeneously over a small volume; in an exact treatment their true distribution on large scales 
must also be taken into account, as discussed in Katz, Balbus & Pacyznski (1986). These authors also identify and describe the 
origin of this logarithmic divergence in the context of microlensing by stars. 

For small r, the sum Ii + h scales as times a term of order unity (due to the logarithmic cutoff) and the phase structure 
function is 



D^{r) ~ AnK^crM-r- = 1.94 X 10"^ (1 + Zifv^ 



M 



IM0/ \l00stai-spc-2/\lpc 



For realistic stellar densities, this implies a very large diffractive scale length. 



'"diif 



2.2 X 10^(1 +Zl)-'v-' 



M 



100 stars pc 



-2 



-1/2 



pc. 



(63) 



(64) 



We note that raiff scales proportional to A^'^^, which implies that the typical scattering angle due to lensing by uniformly distributed 
stars scales as flscat N^^^^. This result has also been derived in the context of microlensing of a homogeneous distribution of 
stars by Katz, Balbus & Pacyznski (1986). These authors point out that if the stellar distribution is not uniform on large scales 
- as must be the case in practice - the large scale stellar distribution also causes macrolensing, which shears the microlensed 
images. This effect is implicitly ignored here because we are only interested in the stochastic fluctuations in the gravitational 
phase delay, and not in large scale phase gradients due to macrolensing. 

The magnitude of r^is is large, which is a reflection of the fact that, in a uniform distribution, each star's contribution to the 
power spectrum is independent of all other stars. The power spectrum of mass density fluctuations scales only linearly with the 
number density of lensing objects, cr, rather than as cr^, when clustering is important (viz. eq. i57\ . Evidently, stellar clustering 
must be important if stars are to make an appreciable contribution to the power spectrum of mass surface density fluctuations in 
a galaxy. 



5.2. Scattering by gas 

Observations of ionized plasma and neutral hydrogen in the local universe suggest that the gas distribution in galaxies follows a 
power law with an index between -3 and -4 from sub-parsec to kiloparsec scales (Armstrong, Rickett & Spangler 1995; Dickey 
et al. 2001; Stanimirovic & Lazarian 2001; Braun 1999). This motivates us to consider a model for the gas whose volume density 
fluctuations follow a power law power spectrum in between minimum and maximum spatial wavenumbers: 

f 0, q< q^nin 

%{q) = I A q-f^, q^in <q< ?max (65) 

\Q, q> q^-ix, 

where the amplitude A is chosen so that the total variance in the volume density is cr^, and we assume q^ni^ <c qmax- For /3 < 3 the 
mass variance is dominated by fluctuations at the inner scale, /o = 1 Iqmax, with A - An{3 - /3) o"p- On the other hand, with 
jS > 3 the mass variance is dominated by fluctuations at the outer scale, Lq = 1/g'min, with A -4n(J3-3) cr^. 
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The two-dimensional power spectrum of surface mass density fluctuations is obtained by projecting the power spectrum of 
volume mass density fluctuations, <l>p, onto the lensing plane. For a gravitational wave propagating along the z-axis one sets the 
argument q, in Op to zero, so the resulting power spectrum of surface mass density fluctuations is ^j,(qx, qy) - AL <^p{qx, qy, q? - 
0), where AL ^ Lq is the thickness scattering medium along the z-axis. (One can think of the quantity <i>p(q_y, qy, q^ - 0) as the 
power in surface mass density fluctuations per unit length.) The phase structure function due to these mass fluctuations is, using 

D^(r)^4nAALK^ q-'-f' [1 - Uqr)] . (66) 

For spectra steeper than /3 > -2 the phase structure function is strongly dominated by mass fluctuations at the outer scale, and 
one has 



DiiAn — q ■ 

2. + p 



l-iFi 



2' ' 2' 4 



(67) 



The phase structure function saturates at the outer scale of the distribution, for r ^ Lq. For small arguments, q^^^r^ - r^/L^ < 1 
and yS > we can expand the hypergeometric function 1F2 to yield 

DJr) = — '-r" X '0^0 [-) , < ^ < 3 ^^g^ 

' P \-Ll, /3>3. 

We identify the quantity s cr^Lg as the total variance in the mass fluctuations encompassed within each cell of size Lo, which 
allows us to express the phase structure function normaUzed to realistic quantities for an L* galaxy. The phase structure function 
takes the form D^{r) - (r/r^iff)^ with, for/? > 3, 



We have normaUzed the difftactive scale length to a line of sight in which the gas fluctuations extend along a distance of 2 kpc, 
while the outer scale of the turbulence, Lq, is normalized to a value suggested by neutral hydrogen observations of galaxies in 
the local group. The mass variance is normalized to a value of only 5 x 10^ M© at the outer scale of the turbulence. A typical 
L, galaxy contains a mass of ~ 5 x 10^ Mq. If most of this gas is distributed inhomogeneously and contained within a typical 
volume ~ 100 kpc^, then the variance of mass fluctuations within each cell of size Lg is = (5 x 10^ M©)^. 

The result for y6 > 3 given in eq. \69\ is the most useful since observations suggest a power law index for gas fluctuations in 
the range 3 < < 4. However, for completeness we include the results for < y6 < 3, for which the phase structure function is 

For instance, the case (5-2 then yields a diffractive scale length of 

-U. i AL V'" i L, 



... = 0.026v-(l..,)-'^^^^^j fej ll^j P'^- '''^ 

As the galaxy would in general possess some inclination to the line of sight, the path length through the gas may vary considerably 
from the value of AL - 2 kpc used in the normalizations here. 

It may be wondered why fluctuations in the gas density of a galaxy, which only account for ~ 5 x 10"^ M© of the mass content 
of an L» galaxy, yield a much greater contribution to the phase structure function relative to the large mass ~ 10" Mq associated 
with stars. The difference arises because of the assumptions underlying the distributions of the two forms of matter. The stellar 
content of the galaxy is assumed to be uniformly distributed and thus explicitly ignores stellar clustering. The assumption of 
uniformity ensures that the power spectrum of mass surface density fluctuations increases only linearly with the stellar density. 
On the other hand, observations show that gas tends to follow a power law distribution of density fluctuations, indicating that the 
gas is clustered in a hierarchy of scales. In this case the power spectrum of surface density fluctuations scales with the square of 
the gas density. We note that the presence of stellar clustering in a population of stars would also cause the power spectrum to 
rise with the square of the stellar density, as demonstrated by equation \51\ . 



5.3. Dark matter fluctuations 

Here we consider the distribution of dark matter associated with the halo of a galaxy. Dark matter is also distributed on larger 
scales (i.e. it is associated with clusters of galaxies), but we concentrate only on dark matter on galactic scales. Irrespective of 
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the dark matter distribution on larger scales, any extragalactic source of gravitational radiation is guaranteed to at least propagate 
through the dark matter distribution of its own galaxy in order to reach an observer near Earth. Both of the foregoing treatments of 
mass fluctuations in § ^5.1l & l5T21 are readily generalized to treat fluctuations in the dark matter, whether it is clumped in star-like 
objects, or clustered on a hierarchy of scales and thus distributed according to a power law. 

We consider the dark matter associated with a galaxy similar to the Milky Way. Here the dark matter must extend to at 
least several kiloparsecs above and below the Galactic plane, and extend out to ~ 30/; ' kpc along the plane. The mass to light 
integrated over the entire Galaxy is ~ 30Mq/Lq; however, this ratio increases sharply near the edge of the luminous disk to 
~ IOOOM0/LO (see, e.g., Sofue & Rubin 2001). 

A plausible model expected on the basis of cold dark matter (CDM) models is a power law spectrum of mass fluctuations 
between some outer scale Lq and some inner scale k with an index (3 ~ 3> for the range of (galaxy-size) scales of interest here 
(e.g. Peacock 1999 and references therein). (The power spectrum is expected to follow a/3 = 3 index for scales below the horizon 
at the epoch of matter-radiation equality ~ 16(/zQ) ' Mpc, where h is the Hubble constant in units of 100 km s 'Mpc '.) 

The model presented in ^5.2l for gas fluctuations is readily adapted to the present situation. As for fluctuations in the gas dis- 
tribution, the two dimensional power spectrum of surface mass density fluctuations is obtained by projecting the power spectrum 
of volume mass density fluctuations onto the lensing plane. One then has ^^.iq.x^qy) - hL<bp{qx,qy,qz - 0), where AL is the 
thickness scattering medium along the direction of propagation, assumed to be along the z-axis. 

The power law scattering for yS ^ 3 is independent of the inner scale, and yields 

We have normalized the diffractive scale length to a total mass ~ 10"Mq distributed inhomogeneously in each cell with radius 
comparable to the outer scale Lom ~ 10 kpc. The outer scale may be larger and encompass a greater mass, but if Mo- increases 
roughly proportional to L^, the diffractive scale remains constant. Gravitational scattering effects in the context of this model are 
expected to be very important, as demonstrated in the following section. 

Another possibility, though somewhat implausible, is that the dark matter is homogeneously distributed in clumps of uniform 
size. In this case the results of 35. H are applicable. However, even for dark matter masses ~ lO'"' Mq, the power associated with 
a uniform distribution of lensing objects is small, and we see from equation \6A\ that the diffractive scale length is many orders 
of magnitude larger than the Fresnel scale. Clearly, gravitational scattering eflects would be unimportant for lensing caused by a 
uniform distribution of dark matter. 



6. Discussion: Implications for the detection of gravitational waves 

We have developed a general theory describing the variability and temporal broadening of gravitationally scattered gravitational 
radiation, and have considered models for the mass fluctuations likely to drive these effects. However, we also need to consider the 
geometry of the lensing situation in order to discuss the implications of gravitational scattering on the detectability of gravitational 
radiation. A large variety of scattering parameters and geometries is expected, depending on the type of gravitational wave source 
under consideration and the astrophysical environment in which it is likely to be encountered. Our aim in this paper is to illustrate 
that scattering effects are important in many instances likely to be encountered in practice; an exhaustive treatment of all the 
circumstances under which gravitational scattering may be important is beyond the scope of the present paper. 

Gravitational radiation always encounters at least two regions of mass fluctuations along its path toward Earth. It must propa- 
gate through the star, gas and the halo of dark matter distributed throughout the source's host galaxy and, to be detected at Earth, it 
must also propagate through the matter distribution associated with our own Galaxy. Propagation through the matter distribution 
of the host galaxy is particularly important for certain sources of gravitational radiation, such as coalescing supermassive black 
hole binaries, as these objects are likely to be located at the centres of galaxies. The radiation is then subject to gravitational 
scattering through a line of sight intersecting the densest environments of the host galaxy. Which of the three constituents of a 
galaxy - gas, stars or dark matter - contributes most to the scattering depends on the distribution of the matter in each of these 
three forms. 

Gravitational radiation may encounter additional mass fluctuations due to any intervening galaxy interposed along the line of 
sight between the source and the detector. Although this eventuality seldom occurs for sources of gravitational radiation in the 
nearby universe, it is worth considering in the context of lensing of powerful sources located at high redshifts, where a line of 
sight is likely to intersect several galaxies and proto-galaxies. For example, the sensitivity of LISA is expected to be sufficient 
to detect coalescence of supermassive binary black holes out to redshifts z ~ 20, and in this case the radiation would likely 
propagate through a large number of intervening systems. 

The geometry of lensing of an extragalactic source due to matter in the Milky Way is similar to lensing by the host galaxy, for 
reasons discussed below. For a source located in the plane of our own Galaxy, the results are quantitatively similar to the lensing 
of an extragalactic source by its host galaxy. 
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6.1. The time scale of wave amplitude decorrelation 

The wave field of a lensed gravitational wave varies due to the movement of fluctuations in the gravitational potential transverse 
to line of sight. The time scale of wave amplitude variations is important, because decoiTelation must occur on a time scale longer 
than v"' for gravitational radiation of frequency v to be detected. The results of ^3. II show that the wave amplitude fluctuations 
decorrelate on a time scale t - rdift /veft , which is the time over which the mean square phase on the lensing plane changes by one 
radian as the phase fluctuations are advected past an observer at speed Veff. This time scale is independent of the distance to the 
mass fluctuations driving the phase fluctuations. 

The velocity of lensing material across the line of sight to the source of gravitational waves is uncertain. The effective speed 
at which the mass fluctuations move across the line of sight depends on the velocities of the lensing medium, the Earth and 
the source. The motion of the source relative to the lensing medium is likely to dominate the effective lensing velocity. Typical 
stellar motions and galactic gas velocities are of the order of one hundred kilometers per second within a galaxy. If such a speed 
is representative of motions of most matter within the galaxy, then it is also representative of the expected speed of any source 
of gravitational waves relative to other constituents of the galaxy. Thus the typical lensing speed is of order v'eff = 100 km s"^ 
Lensing of a gravitational wave source by other material in the source's host galaxy gives rise to wave amplitude decorrelation 
on a time scale 

-.3 / ''diff \ V'eif 



f= 1.8x10^ T days. (73) 

\0.001pcj\l00kms-V ^ 

Temporal decorrelation of the wave amplitude is expected to occur much more rapidly at higher frequencies. The diffractive 
scale is inversely proportional to frequency in the lensing models considered in the previous section. Higher relative velocities 
are expected if one considers lensing by an external galaxy that intersects the line of sight to the source, where the peculiar speed 
of a galaxy transverse to the line of sight is of order 1000 km s"'. 

Scattering-induced temporal decoiTelation of the wave field occurs on a time scale much longer than the period of the radiation 
itself. Thus we conclude that while scattering may alter the amplitude of the lensed radiation, the decorrelation occurs sufficiently 
slowly that the wave field is well correlated from one wave period to the next. 



6.2. Temporal broadening 

Temporal broadening asymmetrically smears the intensity profile of any intrinsically variable gravitational wave source with 
broadening time to - r^/r^^^o). Temporal broadening is important when the ratio rp/rdifr is large, and it depends critically on 
the Fresnel scale, and thus on the lensing geometry. The Fresnel scale depends on the ratio of lensing angular diameter distances 
Deff = DlDls IDs (viz. eq. (|8}) which is well approximated by the source-lens distance when the lensing occurs in the host 
galaxy (i.e. Dls ~ Ds), and by the observer-lens distance when the lensing occurs in our Galaxy. In a Euclidean spacetime (in 
which Ds - Di + Dls the maximum value of Dgff occurs when the lensing plane is midway between the source and observer, 
and DefF = Di/2. However, in our universe the angular diameter distance saturates at z ~ 1 and falls slowly for higher redshifts. 
Thus, for large distances between the source, lensing plane and observer, the slow decline of angular diameter distances beyond 
z ~ 1 implies that the maximum effective angular diameter distance is well approximated by its saturation value, ~ 1 Gpc. 

The thin screen approximation used here is an excellent approximation when treating lensing by an intervening galaxy 
because the thickness of the lensing material is small compared to the effective distance. Deft to the lensing material. However, 
the validity of this approximation is less obvious when treating lensing due to mass fluctuations in the host galaxy of a source, 
where the fluctuations occur over a range of distances that is potentially large relative to Deft. However, the results derived here 
are still applicable provided that one choses an effective distance Deft that is characteristic of the distance to the bulk of the 
lensing material. For instance, in treating the temporal broadening scattering of pulsar radiation due to interstellar scintillation, 
in which the interstellar medium is extended along the entire path from the source to the observer, one takes Deft to be of order 
half the distance to the pulsar. Moreover, a rigorous generalization of the results for an extended medium shows that the scatter 
broadening time may still be written as Deff6'^/2c, where 6c - (^'"diff) ' is the characteristic scattering angle (Lee & Jokipii 1975). 
This argument is applicable to the present situation, as the physics underlying interstellar scintillation and gravitational-induced 
temporal smearing is identical. 

In terms of normalized values typical of scattering in the host galaxy and an intervening galaxy, the Fresnel scale for lensing 
at a distance Deft is 

rp = 0.0028 v-'/2(i+^^)-i/2 /__\ pc^i2v-i/2(i+^j-i/2|__l p^. (74) 

The large value of the Fresnel scale typical for matter located in between the host galaxy and the Milky Way enhances the 
effectiveness of temporal smearing. 

First consider temporal broadening due to lensing caused by the mass distribution of the host galaxy. The diffractive scale es- 
timated for lensing due a uniform distribution of stars is so large that temporal smearing due to scattering by stars is unimportant. 
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The effect of scattering by the gas content of a galaxy depends on the index of its power law distribution. Observations suggest 
an index p - 3 - 4 (e.g. Dickey et al. 2001), so the difFractive scale is of order 0.018 v"'(l + Zl) ' pc for the gas content typical 
of an L» galaxy. Thus lensing due to gas density fluctuations alone causes temporal smearing on time scales of order 4 ms. 

Scattering eff'ects due to the large mass associated with dark matter fluctuations are potentially important. For instance, the 
diffractive scale estimated for lensing by a power law of dark matter fluctuations with an index /3 = 3.2 is r^iff ~ 7.2 x 10""^ v"'(l + 
Zl) With a Fresnel scale rp - 0.0028 v"'^^(l + z/,)"'^^pc, the temporal broadening time is to ~ 2 s. Temporal broadening on 
this time scale would severely decrease the detectability of sources whose intrinsic variations occur on a comparable time scale, 
and render undetectable any intrinsic variations substantially shorter than this scatter broadening time. The temporal broadening 
time is independent of frequency for the models of the matter distribution considered in ^ 

The large value of the Fresnel scale associated with lensing by a galaxy interposed along the line of sight renders the temporal 
intensity variability of most sources of gravitational radiation undetectable. The broadening time for the dark matter difFractive 
scale considered above and the Fresnel scale 1 .2 v"'^^(l pc is a; 4 x 10^ s. The broadening time due to gas fluctuations 

alone is fo ^ 7 x 10^ s. 

The statistical approach adopted in the present treatment of temporal broadening is only valid when a large number of lensing 
lumps of matter - be they stars, gas or dark matter - contribute to the lensing at any instant. This is guaranteed to be the case if 
a large number of objects are contained within a volume of radius rp integrated along the line of sight. However, as the Fresnel 
radius scales as v"'''^, there must be a sufficiently high transition frequency at which the statistical approach becomes invalid. 
Thus, although the diffractive broadening time is independent of frequency, there must be a sufficiently high frequency at which 
the character of the temporal smearing changes. When relatively few objects contribute to the lensing, there is only a small number 
of time delays associated with the small number of ray directions that scatter radiation into an observer's line of sight. An observer 
would then receive a small number of delayed copies of the same intrinsic intensity variations. This transition frequency depends 
strongly on the minimum scale length on which matter exhibits structure. If it is distributed inhomogeneously on scales smaller 
than the Fresnel scale then the present statistical treatment is always valid. The statistical approach is an excellent approximation 
when considering scattering by an intervening galaxy, where the Fresnel scale is particularly large. 

The statistical treatment of lensing in the host galaxy is reasonable if the dark matter is cold. For dark matter particles of rest 
mass energy £„, the spectrum of fluctuations is expected to be cut off at the free-streaming scale (e.g. Padmanabhan 1993), 



which, for cold (E„, > 1 GeVc"^) dark matter particles, is comparable to or below the scales typically probed by gravitational 
scattering in a host galaxy for wave frequencies v ^ 1 Hz. On the other hand, dark matter particles considerably lighter than 
this are expected to exhibit little structure on scales probed by gravitational lensing; a statistical approach would then no longer 
be appropriate, and the character of the scattering would be qualitatively different. In this sense, measurements of temporal 
broadening over a range of frequencies, and thus length scales, provide an exquisitely sensitive means of measuring the nature of 
dark matter. 

7. Conclusions 

The gravitational potentials of mass fluctuations encountered by gravitational radiation as it propagates toward Earth perturb its 
wave front, causing an observer to perceive variability in the wave amplitude and temporal smearing of the signal. We employ a 
statistical approach to relate these effects to the underlying power spectrum of mass density perturbations that drive the scattering. 
This approach is motivated by the fact that gravitational radiation produced by extragalactic sources is likely to encounter the 
gravitational potentials of many objects along its path toward Earth. 

Mass fluctuations along the line of sight to a source of gravitational radiation induce variations in its wave amplitude as it 
propagates toward Earth. The characteristic fluctuation time scale is the time on which the wave amplitude fluctuations decoiTe- 
late. This occurs over the interval over which the root mean square phase difference along the line of sight to the source changes 
by one radian. This time scale is frequency dependent because the phase induced by gravitational perturbations is linearly pro- 
portional to frequency. Wave amplitude variability therefore occurs more rapidly at high frequencies. We consider lensing by gas 
and dark matter, and estimate that decorrelation occurs on a time scale ^ lO^v"' days. However, even when temporal variations 
are too slow to be discerned, the random amplification of the radiation due to focusing and defocusing of the radiation still alters 
the detectability of the signal. 

The most important scattering effect relates to the temporal smearing of scattered radiation, as it has a direct bearing on the 
detectability of any source whose intensity, the square of the wave amplitude, varies with time. The lightcurve of a temporally 
smeared source in a stochastic medium is the intrinsic intensity lightcurve of the source convolved with a function whose width 
is the temporal broadening time scale, to- Temporal smearing arises in a scattering medium because radiation scattered through 
an angle and deflected back towards towards the line of sight takes longer to reach an observer than radiation propagating directly 
along the line of sight. The magnitude of this effect depends on the power spectrum of the mass fluctuations and on the lensing 
geometry. The power spectrum depends on the particular distribution and density of stars, gas and dark matter along the line of 




-4/3 
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sight. Gravitational radiation is always subject to scattering, as it must escape through the mass distribution of its host galaxy and 
propagate through the mass distribution of the Milky Way in order to be detected. 

We consider the consequences of lensing due to a power law distribution of mass fluctuations, due to either gas or dark matter; 
we also consider a simple model due to lensing by a collection of stars. These estimates suggest that temporal smearing due to the 
host galaxy is dominated by the dark matter distribution, and is important for radiation whose intensity varies on a time scale less 
than a second. The temporal smearing time is independent of the observing frequency. Temporal smearing effects scale linearly 
with the effective distance to the scattering material, Defr = £>l£>ls IDs , so the smearing time due to lensing by an external galaxy 
is exceptionally large. For lensing at cosmological distances, Deff ~ 1 Gpc, the temporal smearing time is estimated to be of order 
to ~ 10^ s. As the distribution of dark matter on small scales relevant to temporal broadening is ill-constrained, we caution 
that the foregoing estimates are highly uncertain. When dark matter is absent, even the gas content of a galaxy can appreciably 
broaden the intrinsic intensity variations of a source. The broadening time associated with the gas content of the host galaxy is 
of order milliseconds, while the time associated with scattering due an intervening system is of order 10^ s. 

The circumstances in which temporal smearing poses a serious limitation on the detectability of intrinsic source variations 
depends on the distribution of scattering material. The description of temporal smearing presented here is valid when one is 
justified in treating the phase fluctuations using a statistical approach. However, it fails when only a few inhomogeneities con- 
tribute to the scattering at any one instant, at which point the observed intensity lightcurve becomes the sum of a small number of 
delayed copies of the intrinsic intensity lightcurve. Thus its validity varies according to the line of sight and the scale on which 
matter distributed along it is inhomogeneous. Clearly, the scale on which dark matter is inhomogeneous is critically important. If 
it is inhomogeneous on scales smaller than the Fresnel scale, as is expected for CDM fluctuations, then the statistical approach 
adopted for temporal broadening here is well justified. Thus, the scale on which dark matter is distributed is of critical importance 
in determining the detectability of gravitational radiation. Conversely, temporal smearing measurements of gravitational radiation 
place strong constraints on the distribution and hence the nature of dark matter along any line of sight through which a source of 
gravitational radiation is detected. 

Acknowledgements. The author thanks Don Melrose, Mark Walker and Ole MoUer for their encouragement and comments, and Marco Spaans 
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Appendix A: Derivation of mutual coherence function V{ki , A;2; r = 0) 

Here we calculate the mutual coherence function, {4>(k+Ak/2, z, r)^*{k-Ak/2.z, r)), following on from equation i39\ . We perform 
the average over phase fluctuations assuming that iff is a normally distributed random variable with zero mean (the mean phase 
delay is incorporated into the term involving the group velocity in equation (I36» : 
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With the change of variables 
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which has a Jacobian 1 + 0(Ak /k ) which we approximate as unity, the mutual coherence becomes 
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We make the further change of variables r' = r + sAk/2k and s' = s, which again has a Jacobian 1 + 0{Ak^/k^), to write 
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A further change of variable q = s'/'"p yields the final result. 
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This expression is equivalent to that obtained by Lee & Jokipii (1975) for interstellar scattering through a medium of plasma 
inhomogeneities. 



Appendix B: The effective velocity of scattering material across the line of sight 

We derive an expression of the velocity of scattering material across the line of sight, taking into account motion in the source, 
lensing medium, and the Earth. Consider an object with angular diameter distance Ds from Earth moving at a transverse velocity 
Vsic measured relative to some frame. Between two time intervals t - Q and t - the object traverses a displacement Vsicfo, and 
the Earth's change in position is VEarthfo, where VEarth is measured relative to the same co-ordinate frame. From Fig.|Blwe see that 
the line of sight moves a distance As = VEaithfo + (Vsrc - '^Eaiih)t(i{DLl Ds) across the scattering screen. The effective velocity of 
the line of sight across the scattering screen is thus As/Jq: 



Veff = VEa..h|l-f^)+V.J^ 



(B.l) 



In this time interval the screen also moves an additional distance Vscjeenfo. so the eff'ective velocity at which the line of sight to the 
object crosses a point on the screen is 



Veff ~ ^screen ^Eaith \^ I ^^rc I f-. 

\ L>s I \Ul 

The observed scintillation time scale is the time taken by the line of sight to traverse a length scale rjiff on the screen: 

'"diff 



V'eff 



(B.2) 
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This is the time scale of the scintillations as observed at Earth. 

There is a subtle distinction between the velocity at which the source-observer line of sight moves across the lensing plane 
and the actual velocity at which an observer moves through the lensing pattern, as it is projected at Earth. The distinction arises 
because the wavefronts emitted from the source are spherical, and an observer located near the lensing plane would not measure 
the same lensing pattern scale length as an observer on Earth. A wavefront "imprinted" with a pattern of scale length rdiff.screen on 



Jean-Pierre Macquart: Scattering of Gravitational Radiation 



19 



Vsrc to 




VE to 

Fig.B.l. The geometry showing the distance traversed across the lensing plane in a duration fo. 



the scattering screen has a pattern scale of length rdiff,Eaiih = '"ditf.scieen^'s IDis upon reaching Earth. Thus, for an observer who 
measures rjiff ,Eaiih directly from the scale size of the scintillation pattern passing across the Earth, the real speed of the scintillation 
pattern at Earth is 

^diff ,Eaith ,t-, , . 

Vreal = • (B.4) 

^scint 

Thus one has 

(^)'-diff,scieen Ds Ds - 
Vreal - nirsc --n ~ '^'screen — VEaith - 77 Vsrc- (.D.O) 

.; L)ls Uls Uls 

This correction is only a concern if we were to measure the scale size of the lensing pattern directly with an interferometer on 
Earth. The expressions derived above reduce to those derived in Appendix C of Gupta, Rickett & Lyne (1994) for scattering in 
Euclidean space. 



